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伪 抛 物 型 积分 微分 方程 的 混合 有 限 元 误差 估计 ， 


车 海 涛 
(山东 省 潍坊 学 院 数学 与 信息 科学 学 院 ， 潍 坊 261061) 
摘 要 : 基于 Raviart-Thomas 空间 ， 本 文 对 伪 抛 物 型 积分 微分 方程 初 边 值 问 题 提 出 了 混合 有 限 元 方法 。 
与 通常 的 有 限 元 方法 相 比 ， 该 方法 可 以 同时 高 精度 逼近 未 知 函 数 及 未 知 函 数 的 梯度 。 通 过 引入 广 
义 混合 椭圆 投影 ， 证 明了 其 存在 唯一 性 ， 并 得 到 了 其 一 系列 性 质 ， 利 用 其 性 质 给 出 了 平方 模范 数 
下 的 最 优 误差 估计 ; 利用 广义 混合 椭圆 投影 和 正则 Green 函数 得 到 了 最 大 模范 数 下 的 拟 最 优 误差 


估计 。 
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1 引言 


SIER? PRE Lipschitz 连续 边 界 6Q 的 有 界 区 域 ,对 固定 的 满足 0 < 了 < ce 的 T， 我 们 
考虑 如 下 的 伪 抛 物 型 积分 微分 方程 初 边 值 问题 
Ut = div(a 7 w +b 了 以 十 fy be Vv udr) +f, (2,t)EQ~x [0,7], 
u(x,t) =0, (a,t) € AQ x [0,7], (1) 
u(x, 0) = uo(z), LEQ 
WA TATE, PR BMa(z,t), bi (x,t), bas, t 以 及 它们 的 导数 光滑 有 界 且 满足 0 < co < 
a(z,t) < co WERI FERH f € C1(0,T; ZL2(Q))， uo € 五 *(Q)，(s 是 一 固定 常数 )， 问 
题 (1) 都 是 唯一 可 解 的 。 


令 


t 
p= —a Vu: — bı Vu 一 J boVudr, 
0 
并 设 
a(z,t)=a™"(z,t), (x,t) = a(z,t)bi(z, t), C(x, t, T) = a(x, t)be(z, t, T), 


B(x, t) = —Vb(z,t), yl(z,t, T) = —Ve(2,t,7), 
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则 问题 (1) 可 写成 如 下 的 混合 一 阶 系统 
ut + divp = f, (z,t) EQ x (0,T], 


ap + Vu + V(bu) + Bu + IN V(cu)dr + f yudr =0, (2,t)€ 0x (0,TI, (2) 


u=0, (x,t) € AN x (0,T], 
u(0) = uo, tE Q, 
TERE uls, tloa =0, RIRE ulz, t)lan = 0。 
定义 空间 


V=H(div,0)= {v € (L?(Q))"; divv € D7(0)}, W =N), 
其 中 V BEA (vi = jivi + |ldivv|]?. W (2) 的 弱 形 式 为 : 求解 {p,w} : [0, 相 一 V x WW， 使 得 


(ut, w) + (divp, w) = (f,w), Vwew, 0<t<T, 

(ap, v) — (ue + bu + fy cudz, divv) 加 
+ (Bu + fj yudr,v) = 0， vveV, 0<#t<T, 

(u(0), w) = (uo, w), VweW. 


为 了 定义 {p, 好 的 一 个 适当 的 有 限 元 逼近 过 程 ， 考 虑 将 9 剖 分 成 三 角形 单元 的 拟 一 致 部 
4} Th (其 单元 直径 不 大 于 h(0 < h < 1)) 相 联 系 的 有 限 维 子 空间 VW x Wro HV, x Wh C 
V x W, Th 的 边界 单元 允许 有 一 条 边 是 曲 的 。 我 们 选择 Vi x Wi, 为 Raviart-Thomas ZP, 
其 指标 大 > 0。 引 进 世 投影 ph : W 一 Wi M Raviart-Thomas#t# Ip : HHQ)? 一 Vh, BRA 
交换 性 质 div o II, = pp o div: H! (Q)? — Whe 

方程 (3) 的 关于 时 间 连 续 的 混合 有 限 元 逼近 为 : 求 {pi,un} 一 Vax Wh， 使 得 


(unt, w) + (divpn, w) = (f,w), VweWh, O<t<T, 

(apn, V) — (une + bur + i cundr, divv) 7 
+ (Bun + fo Yundr,v) =0, VvEVn, 0<t<T, 

(un(0), w) = (uo, w), Vw € Wh, 


且 (4) 式 的 解 存 在 唯一 。 
本 文 余下 部 分 安排 如 下 : 在 第 2 节 将 讨论 与 (3) 式 相 联系 的 广义 混合 椭圆 投影 的 性 质 ， 在 
第 3 节 将 展示 本 文 的 主要 结果 。 
2 广义 混合 椭圆 投影 
引入 与 (3) 式 有 关 的 广义 混合 椭圆 投影 {Pr tn}: [0,T] 一 Vax Whp， 使 得 
(qiv(P — Pah w) =0, Vwew,, O<t<T, 
(a(p — Pa), v) 一 ((ue — tnt) + b(u — Un) + Rk c(u — up)dr, divv) 
+ (Blu — ŭn) + Ralu — ŭn )dr, v) =0, Vvev,, 0<t<T, 
(uo — ŭa (0), w) =0, Vwe Wh. 
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E (5) 的 解 存 在 唯一 。 
& pi =p- Ph, U1 =phu — Ùh, U2 =u- phus M (5) KARA 


(divpı, w) =0, Vwewn, 0<t<T, 
(apı, v) — (ure + blui + u2) + IH c(u1 + ug)dr, divv) 

+ (Blur + u2) + fo yur +u2)dr, v) =0, VveVn, 0<t<T, 
(ui1(0),w) =0, Vw € wh. 


引 理 2.1 车 {pi,wi} 满 足 关系 式 (6)， 且 假定 是 2- 正 则 (2- 正 则 性 的 定义 见 文献 凤 )， 则 
对 所 有 的 0 <t< 人 下， 存在 不 依赖 于 h 和 t 的 常数 C > 0， 使 得 


WESC 
leas ell < C (hllpill + hlluall + lasl- + he divp | 
+ f (hipa + 42% ldivpall + hull + lluall-1)47), 
其 中 当 k 二 0 时 ，6ko =1; SKD1IN, ôk = 0。 


引 理 2.2 $ {p u} M {pp, dn} DHE (3) 和 (5) KH. BR {p u} HAI ANKE, 
则 对 0 <t < 人 下， 存在 不 依赖 于 及 和 t 的 常数 C > 0， 使 得 天 之 0， 且 


lidiv(p — Ba) < Chrlpllrrs, 0<r<k+1, 


t 
lonu — ŭa l| < Chto f ell- + lull-)dr, 1<r<k+1, 
0 


t 
pre — tne CRS {Hull + elle + f (pl + lular}, Vr Sk +1 
0 


Ip — Ball < caflph + lula + f (pth + Huldar}, 
Ju Gl < Ohf fal + f(s + lular}; 


t 
lue — Gell < Ch{ falls + lls + uelli + f (pla + lear}, 
0 
2) 4k>1,2<r<k+1h, 


lp ~ Ball < CA"{Ipll + lle + f Opl- + Heda} 
[u tal < CH" {llul + f (pl-i + lule—a)er} 


ue — Binal < CA {ll + ema + lela + f (ella + Hls)2r}- 
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下 面 考虑 {Pa tn} 的 Zee(Zece) 误差 估计 。 引 入 如 下 定义 的 正则 Green 函数 外 


G,+VA,=0, TER, aG + VrA2 = 68, TER, 
divG1 = 8È, TER, 和 divG = 0, TER, 
A1 = 0, ZEO0 A2 = 0, z EQ, 


HH Of, 小 是 z € OM) Dirac RR, HWE 
(w, ô?) = w(z), VweWn,  (v,d2)=v(z), Vv eu. 
对 于 z ARASH 
leello,c < 2I(w, 69) 1vllow < 2I(v, ô). 


令 {GT AT € Va x Wn Æ (Gi, Ar} MRA APR TCE, {G}, Ak} © Va x Wi Æ {G2, ro} 的 
混合 有 限 元 逼近 。 由 文献 由 和 [5]， 我 们 得 到 


AZ + GE < Clin pl 二  ||G8llo1 < Clima. 
引 理 2.3 在 引 理 2.1 的 条 件 下 ， 对 0 < 上 < 了 和 0 < h <1, 则 存在 不 依赖 于 疡 和 + 的 党 
数 C > 0， 使 得 
1) 4k-On, 


t 
lu — rllojoo < Ch| Inh]? {ul + f (lulla, + lliplli Jar}, 
i t 
lte 一 各 om < Chl In Mè {leo + lulla + hells + f (lls. + pla Jar}, 
0 


t 
Ip = Balloo < Chl nl {Ilo + plio + f (luli + ell Jar} 


2) 4k>1,2<r<k+1m, 


t 
[u Talos < CH" {al + f (lul + pli )ar)}, 
” t 
fu: = Un,tllo,co < Ch" { luclln,oe 十 lull, 十 llplir + i ( llull- ag ilipli- Jdr}, 
t 
IP = Baloce < CA mhl {Hulse + Helse + f° (lle + Ile ar}. 


3 ”连续 时 间 的 混合 有 限 元 逼近 


首先 利用 广义 混合 椭圆 投影 {54,i} 得 到 {pn, un} 的 最 优 L? 误 差 估计 ， 然 后 利用 投 
影 {Bn, Unt 和 正则 Green 函数 得 到 {Ph, un} 的 拟 最 优 L” 误差 估计 。 

定理 3.1 令 {p,zj Hl {pr un} 分 别 是 弱 形式 (3) 和 离散 形式 (4) 的 解 ， 假 设 p, u, 充分 光 
消 ，0 是 2 正则 ， 则 对 0 < t < T， 存 在 不 依赖 于 及 和 t 的 常数 C > 0， 使 得 
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1) 当 k=0 时 ， 
lu un] < Cn lial + f (vals + lips + afar}, 
llus — un,ell < Chf [lull + |lwella + lliplli +f (walla + liplls + lulyar}, 
Ip = pall < Ch{ lal + uelli + bls + f (uelh + pls + lallar}, 
0 


t 
ldivp = pal] < Chf llulla + eel + ela + f (bull + pll + falar} 


2) 4k>1,2<r<k+1M, 
t 
Juan < Oh" {halle + f (alle + liple + lhale—a)4r 
t 
ue — uneli < CH" $ ala + lel + ella + f (lel + Hella + lle)ar}, 
0 
t 
Ip — Pall < CA {alles + [fuels + pl-i + f (uelle + plea + [lull--a)er}, 
0 


t 
laiv(p — pa) < CH” {hls + luelh + linea + f (ueli + lle + ella}. 


定理 3.2 在 定理 3.1 的 条 件 下 ， 对 0 < t+ < 工 和 0 < h < 4, MEEDRKIRF hth 
数 C > 0， 使 得 


t 
le wall < CII hl {uo + f Cease + ell + lll)ar), 
t 
ue — thello,co < Chl In A? {lleo + Iftulla,oo + lplli + I (lelli + llulla, + plli)ar)}, 


t 
lip — pr llo,so < Ch In h]} {ulli,co + IIPli,co + f (uella + llulls.oo + Ilpll)a7 }- 


2) 当 大 > 1 2<r< 大 十 1 时 ， 得 到 
i t 
lu — uno < CR] im A {ln + f Celis + lelle + Iplar} 
0 
4 t 
l — unloeo < CHin hl {else + Hell + elle + f, (lel, + Hell + plar}, 


i 
3 
IP = Palace < Ch" IIn h3 {lulose + plno + 人 del + lle + Iplar}. 
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Error Estimates for Mixed Finite Element Methods for 
Pseudo-parabolic Intergo-differential Equations 


CHE Hai-tao 


(College of Mathematics and Information Science, Weifang University, Weifang 261061) 


Abstract: In this paper, a mixed finite element method is proposed to solve the initial-boundary 
value problem of pseudo-parabolic integro-differential equations. Compared with the usual finite el- 
ement method, the unknown scalar and the adjoint vector function are approximated optimally and 
simultaneously. By introducing the projection of generalized mixed element, the existence and unique- 
ness of projection and some important properties of the projection are proved. By using its properties, 
the optimal order error estimates in square norm are derived; By using the projection of generalized 
mixed element and the regular Green function, quasi-optimal order error estimates in maximal norm 
are finally obtained . 

Keywords: pseudo-parabolic intergo-differential equation; mixed finite element method; error esti- 
mate 


